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Abstract. In this work we study how quantum fluctuations modify the quantum 
evolution of an initially classical field theory. We consider a scalar <fi 4 theory coupled 
to an external source as a toy model for the Color Glass Condensate description 
of the early time dynamics of heavy-ion collisions. We demonstrate that quantum 
fluctuations considerably modify the time evolution driving the system to evolve in 
accordance with ideal hydrodynamics. We attempt to understand the mechanism 
behind this relaxation to ideal hydrodynamics by using modified initial spectra and 
studying the particle content of the theory. 
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1. Introduction 

One of the outstanding theoretical problems in heavy-ion physics is a first principles 
understanding of the isotropization and thermalization of the matter produced in 
collision. The fact that the system is nearly thermal and isotropic at early times has 
been deduced from hydrodynamic model fits [U El El HJ El E] to the measured spectra and 
elliptic flow [HI El HO]- Estimates of the relaxation time range from r re i ax ~ 0.5-2 fm 
[3] which is hard to accommodate within a simple picture of interacting quasi-particles. 

However, a quasi-particle description is not essential to thermalization and in this 
work we will demonstrate that an initially strong classical field undergoing quantum 
evolution may evolve in accordance with ideal hydrodynamics. We will show that the 
presence of secular divergences (modes whose occupation number grows with time) 
become semi-classical on relatively short time scales and must be resumed to all orders 
in a standard perturbative expansion. 

While this works focuses on scalar 4 theory it is suggested that similar mechanisms 
may be at work within the framework of the Color Glass Condensate (CGC) description 
[TT1 [T2l IT3l [T4l [TBI IT6l IT7I ITS] of high energy nuclei. Work in applying the techniques 
shown here to the case of classical Yang-Mills is currently in progress (THJ. The goal 
of this work is to understand the role secular divergences play in modifying the time 
evolution of the classical field. We will demonstrate that the quantum evolution of a 
scalar field evolves in accordance with ideal hydrodynamics. Finally, we will speculate 
on the mechanism behind this relaxation by looking at modified spectra and the time 
evolution of the number density. 

This work is largely based off the first paper on this topic [20] which we refer 
the reader to for more details on the resummation scheme. It is worth pointing out 
that the resummation of secular diverges is qualitatively similar to the resummation 
of leading logarithms (g 2 Xi^)) n of the incoming partons' momentum fractions 
required for the computation of inclusive quantities at leading log order [211 1221 123] . 
These results have proven to be valuable in a quantitative understanding of the near 
side angular correlations observed in nucleus-nucleus J2U [25] and proton-proton [261 121] 
collisions. We would like to point out that a considerable amount of work has been 
done for theories similar to 4 in the context of reheating after inflation [281 EHl 30J. 
In addition, considerable progress has been made on the thermalization problem in the 
context of Heavy-Ion collisions through the use of N-particle irreducible effective actions 
[321 [331 EH E5l [361 [37] . 

2. The Model 



The CGC inspired scalar theory model has the Lagrangian 
£ = \ (d,<f>) - V{(f>) + J<p 



(1) 
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where the interaction potential is 

V(d>) = f^ 4 (2) 

and J is an external source which mimics the large x color charges of the incoming 
nuclei. Since the external current vanishes after the collision takes place we take our 
source to be nonvanishing for x° < only, 

J(x) ~0(-x )^. (3) 

The role of the external source is to initialize a classical field (having occupation 
number ~ l/g 2 ) which evolves solely via their self-interactions at x° > 0. This source 
term also brings an external scale into the problem. Since we assume the source is 
turned on adiabatically from x° — > — oo we can consider the evolution of the classical 
field as an initial value problem at x° = with = 0) ~ Q/g and = 0) = 0. 

Throughout this work we will use the same model parameters. In order to avoid 
confusion we now state these parameters once and for all. For both the homogeneous and 
non-homogeneous systems we take 0o = 12, <fio = and g = 0.5 (a very weak coupling 
considering the factor of 4! in front of the potential). For the case of the 3D simulations 
we employ a 12 3 lattice with a volume of 12 3 . This lattice size will have a momentum 
cutoff of k max ~ 5.44. For these model parameters the resonance mode exists between 
3 ^ k res < 3.22 at t = 0. The zero-mode has an effective mass m 2 = (g4>o) 2 /2 = 18 and 
a period of oscillator of T w 3. 

3. Homogeneous System & Homogeneous Fluctuations 

In this section we consider a classical background field that is homogeneous in all space 
and undergoes quantum evolution with space-independent fluctuations (i.e. zero-mode 
fluctuations). While highly unrealistic this simple toy model will allow us to see how the 
mechanism of phase decoherence leads to ideal hydrodynamic evolution. The Lagrangian 
for a uniform non-expanding scalar theory is 

£ = \tf - VM , (4) 

where = d<p/dt. The classical evolution can be found in closed form. Since the energy 
T-L = 2 /2 + V((p) remains constant throughout the evolution we can write 

^ 2 = E - V(<j>) , (5) 

where Eq is the initial energy of the system which is determined by the initial condition 
of our classical field 

Eo = \& + ^ti ■ (6) 

where O — 4>{t = *o)- Equation [5] can be integrated to obtain 

1 /■*(*) dip 

t-to = ~j= 1 ■ (7) 
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At this point it will be useful to introduce some notation. Let us define 

e 2 = <?74! , (8) 

and make the change of variables yfeip = —E^ 4 cos6. We are also free to set t = for 
the non-expanding case and we find 

t =-L m r , a ( 9) 

where 



i (Vt<Kt)\ 



e(t) = cos 



V El" ) 



e ° =cos "(iK) (io) 

The above integral equation can be solved for <j)(t) in terms of the Jacobi elliptic function 
of the first kind having elliptic modulus 1/2. 

<P(t) = -^cn 1/2 hyfeE^t - F 1/2 (0 O )1 (11) 
V e 

where Fi/ 2 (0 o ) is the incomplete elliptic integral of the first kind of modulus 1/2. The 
above result is periodic with period 

T = ^T 7 I K ( 1 /2) (12) 

where if (1/2) « 1.85407 is the complete elliptic integral of the first kind. Notice that 
the period of oscillations depends on the initial conditions through E . The fact that 
the period of oscillation depends on the initial condition is a signature of non-linear 
evolution and is crucial for phase decoherence. 

It is worth noting that to a very good approximation (within about 15%) the above 
expression for can be approximated by 



E 1/A 

<f)(t) « — j=r COS 



(13) 



(14) 



where £ is a phase set by the initial conditions 
3.1. Stress-energy tensor 

With an analytic expression for <fi(i) available we can now find analytic expressions for 
the stress energy tensor as well. For the homogeneous non-expanding system there are 
two independent components of the stress-energy tensor 

T 00 = \tf + V(<f>) 

T ij = 6* Q</> 2 - V{<t>)) (15) 
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with all other components vanishing. Using the expressions derived in the previous 
section we find 



T 00 = E 
T ij = E 

~E 



1 - 2cn? /2 (2^o /4 i - F V2 ( 
l-2cos 4 [—(*-£) 



(16) 



As an example, in Fig. [T]we plot T 00 and T 11 as a function of time. We also show 
the good agreement between the true solution and approximate form of T 11 in this 
figure. Clearly, this LO result does not have a well defined equation of state. 

3.2. Spectrum of fluctuations 

We now want to superimpose quantum fluctuations on top of our classical background 
field. In this section we will consider the following toy model for the spectrum of 
fluctuations, 

F(a,a) = 5(a)- ? i=exp(-^ , (17) 

where a characterizes the variance of the zero-mode fluctuations. For this toy model we 
will treat a as a free parameter. It will be computed from first principles later on. We 
should stress that this is a highly unrealistic model since we are ignoring any quantum 
fluctuation which are non-homogeneous in space. 

The expectation value of an inclusive operator (such as the stress energy tensor) is 
defined as 

/+0O 
da da F(a,a)Oi,o((fio + ct,<fio + a) (18) 
-oo 

where Clo(0o + a , 0o + a) is the operator of interest computed at leading order with 
initial conditions shifted by a and a. For this particular choice of fluctuations and using 
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Figure 2. T 00 /3 and T 11 for a uniform non-expanding cj) 4 theory after averaging over 
a Gaussian distribution of fluctuations in </> having a = 0.4. The dashed curve shows 
the envelope as given in the text. 



our approximate solutions for </>(t) found in the previous section, the integrals over a 
and a can be done analytically. The result is shown in fig. [2] for o = 0.4. The analytic 
expression is not too enlightening. It essentially consists of a number of terms having 
oscillations at different frequencies which die off exponentially at different rates. But it 
is instructive to pull out the one term which dies off slowest. Its envelope is given by 

~ e - 2c2 ° w , (19) 

where we have defined the constant c = — z m « 0.3459. We can now identify a 

K(l/2)Vi\ 

relaxation time 



1 2 

Trclax = -7= « • (20) 

y/Zcga go- 
While this is a very unrealistic model it is nice that the above result could be derived 
analytically. In the example shown in Fig. [2] the fluctuations are completely absent after 
2 x r rc i ax . In other words, by t ~ 4/ (ga) ~ 20 the system has a well defined equation of 
state (e = 3p) and evolves in accordance with ideal hydrodynamics. 

The mechanism behind the relaxation of the pressure is quantum decoherence which 
we now explain. Each initial condition in our ensemble average is shifted by a random 
Gaussian variable (0o — > <Po + a ) and this corresponds to a shift in the initial energy £o 
of the system. The time evolution of each individual system is periodic with a slightly 
different period of oscillation as given by Eq. [12j When performing the ensemble average 
the differing periods of each system results in a phase decoherence forcing the pressure 
to relax to its equilibrium value. Let us stress that this will not occur in a </> 2 theory. 
In this case the period of oscillation will not depend on the initial condition. 

4. Non— Homogeneous Fluctuations 

In the previous section we showed how a homogeneous system undergoing zero-mode 
fluctuations relaxes to a system evolving according to ideal hydrodynamics. While the 
previous case is of pedagogical interest since it shows simply how the decoherence of the 
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Figure 3. Evolution of a linear perturbation on top of the background field. Left: 
Linear growth of the k = mode. Middle: Exponential growth of a resonant mode 
k r es = 0.5165<7</>o. Right: Typical behavior of a stable perturbation k = 2k res . 



quantum field leads to the relaxation of the pressure it is highly unrealistic in that it 
does not include space dependent fluctuations. 

We now consider the same model in three dimensions including the space dependent 
fluctuations as predicted from quantum field theory. In order to motivate the need for 
the resummation we first discuss the case of linearized perturbations. 



4-1. Linear perturbations 

In this section we now consider how a linearized perturbation evolves on top of the 
homogeneous background field. We decompose the background field into a homogeneous 
part 0k=o and a small field perturbation a{x). The equation of motion for the Fourier 
transform of our field perturbation a(x) is 

a ±k + [k 2 + ^"(0 k=o )]a ±k = O. (21) 

In the above expression </> k= o is the zero-mode solution given by Eq. [Til We now 
numerically solve Eq. [21] in order to investigate how linear perturbations evolve on top 
of the background field. In Fig. [3] we show how the amplitude of three k modes evolve 
when given an initial amplitude of a k (t = 0) = 0.1. The first plot shows the zero mode 
whose amplitude grows linearly with time. The second mode is taken from within the 
resonance band and it clearly grows exponentially with time. The third mode shows 
the typical behavior of a high momentum mode (here shown for k = 2k res ). While the 
high momentum modes can be treated perturbatively as their amplitude does not grow 
with time the lower k modes lead to secular divergences. Clearly, at times when gt 
(for modes outside the resonance band) or ge^ (for resonance modes) become of 0(1) 
a resummation becomes necessary. 



5. Results from the full fluctuation spectrum 

5.1. Initial Condition 

As is clear from the previous discussion quantum fluctuations on top of the homogeneous 
background field will play an important role in the resulting dynamics. The spectrum 
of these quantum fluctuations are derived from first principles. In this case the classical 
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field and its conjugate momentum n = dC/dcf) = <p are promoted to quantum operators 
(f) and 7T obeying the equal time commutation relations 

0(x),#(y)] = i5 3 (x-y) (22) 
The field operators can be rewritten in terms of creation and annihilation operators 
i, s If d 3 k 



(2vr) 3 / 2 7 



+ a k e" 



-ik u x^ 



(23) 



obeying 

a k ,a{J = 5 3 (k-p). (24) 

Using the above mode decomposition one can easily show that the two-point correlation 
function in a homogeneous background field takes the form 

1 r d 3 k 



(0(x)0(y)) 
(7r(x)7r(y)) 



(2tt) 3 
1 

(2tt) 3 



2o> k 
d 3 k 



3 ik-(x-y) 



-w k e 



ik-(x-y) 



(25) 



(26) 



where w k = k 2 + m 
random field having power spectrum 



2 . Our semi-classical simulation will therefore consists of a Gaussian 



2(27r) 3 w k (27) 
superimposed on top of the homogeneous background field. The power spectrum as 
written above is UV divergent and this is regulated by the lattice spacing. If we impose 
a momentum cutoff A the energy density will contain terms that behave parametrically 
as Q 4 /g 2 ,Q 2 A 2 and A 4 . The A 4 is a pure vacuum contribution and can be computed 
by performing simulations with the source J turned off which can then be subtracted 
from the corresponding result. The Q 2 A 2 terms in not renormalizable in the usual sense 
since it mixes diagrams having an arbitrarily high number of loops. In practice, what is 
done, is to choose a cutoff that is sufficiently large in order to encompass the relevant 
physics (A > m) but small enough to keep the cutoff-dependent terms negligible with 
respect to the classical contribution (A Q/ '^/g). 



5.2. Results 

Fig. HJshows the pressure and energy density (e/3) as a function of time with an ensemble 
average of 1000 simulations. The main conclusion of this paper is that the ensemble 
averaged pressure relaxes towards e/3 and therefore has a well-defined equation of state 
and evolves in accordance with ideal hydrodynamics. We should stress that even though 
there exists a well-defined {i.e. time-independent) equation of state the system is not 
necessarily in thermal equilibrium as well will show. This rapid establishment of an 
equation of state known as prethermalization has been studied in the context of a linear 
cx-model using the 2PI effective action in [3H]. This work found similar conclusions; 
there can be the rapid establishment of an equation of state via phase decoherence 
regardless if scattering processes can thermalize the system. 
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It is apparent from figure H] that the time evolution evolves in two stages. First, in 
the window < t < 50 the amplitude of the pressure oscillations decrease very quickly 
to moderate values. Then from a time t ~ 50 and onwards there is a slight rebound and 
a gradual approach to complete relaxation. 

In order to try to interpret this result and understand the role of different excitations 
we perform additional calculations using a modified spectrum of fluctuations. Even 
though these modified spectra will result in the incorrect quantum expectation values 
the results may serve useful in understanding the role of different fluctuations. In 
Fig. [5] we show the resulting pressure after an ensemble average of 250 configurations 
for various initial spectra which we now discuss. In one case (upper left figure) we omit 
quantum fluctuations of the zero mode. In a second case (upper right figure) we omit 
any initial fluctuation within the resonance band. The lower two figures show spectra 
which omit any initial fluctuation having k < 4.4 (lower left) and k > 2 (lower right). 
Of course, once the time evolution begins, there is nothing stopping self interactions 
from causing excitations to scatter into these initially unoccupied modes. Not including 
the resonance modes in the initial spectrum of fluctuations (as done in the top right of 
Fig. |5]) is different from the analysis of [20] where the lattice cutoff was chosen to be 
below the resonance band. In the latter case, the resonance modes can never become 
occupied, which was found to significantly modify the evolution of the pressure. 

There is a lot one can infer from Fig. [5j The upper two figures comprise results 
missing a small fraction of the initial spectrum. In the first case we neglected initial 
fluctuations of the zero mode while in the second case we are neglecting the very few 



Role of quantum fluctuations in a system with strong fields 



10 




Figure 5. Evolution of the ensemble averaged pressure for various incorrect spectrum 
of fluctuations. Each case consists of an ensemble average over 250 configurations. 
Top Left: Omitting the zero mode. Top Right: Omitting modes within the resonance 
band. Bottom Left: Omitting modes with k < 4.4. Bottom Right: Omitting modes 
with k > 2. 



modes that sit within the narrow resonance band. In both of these cases the results 
are qualitatively similar to the results using the full spectrum shown in Fig. HI In the 
case of the zero-mode it is not surprising that its absence doesn't affect the result. 
One can estimate the relaxation time from the zero-mode alone based on the previous 
section where we showed in Eq. [20] that the relaxation time in the homogeneous case is 
inversely proportional to the standard deviation of the Gaussian fluctuations as given 
by the power spectrum. In this case 

r relax « 2 « 180 (28) 

gy/Vt(k = 0) 

for g = 0.5. Clearly, this is a much longer relaxation time then observed in the full 3D 
simulation. In the case of the absent resonance modes, one can see by looking at the 
occupation numbers, that the occupied modes are very quick to scatter and perturb the 
initially unoccupied resonance band. 

What is more interesting is if we neglect a large portion of the initial spectra. 
The lower left plot shows the ensemble averaged pressure with a spectra including the 
intermediate momentum modes (basically we include modes higher than the resonance 
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band). The evolution from < t < 50 is remarkably similar to the result using the 
full spectrum. In the lower right plot we have used a spectrum consisting of only low 
momentum modes. In this case we no longer have the rapid relaxation at t ~ 50 but 
instead have a gradual relaxation that extends to t ~ 150. Based on this analysis 
we can understand the two-stage relaxation observed when using the full spectrum. 
The first relaxation in the period < t < 50 is clearly controlled in some manner by 
quantum fluctuations above the resonance band. While the second (more gradual) stage 
of relaxation, taking place for t > 50 is controlled by modes below the resonance band. 
The time scale for relaxation due to the low momentum modes is on the order of that 
estimated in Eq. [28] for the zero-mode. Of course this interpretation is only qualitative. 
Self-interactions immediately cause modes which are initially unoccupied to become 
occupied and the result becomes a complex interplay between many modes which cannot 
be understood simply by studying the linear evolution of individual quanta. 

It is interesting to note that the time scale for prethermalization in the three 
dimensional simulation has the same order of magnitude (it relaxes about 3-4 times 
faster) as the crude estimate of equation [28] found for the homogeneous case. If we take 
the estimate from eq. [28] seriously we see a faster relaxation with increasing coupling 
constant g. The power spectrum entering into eq. [28] is determined by the quantum 
field theory in question. Modes with lower k will have the largest fluctuations. In our 
case V^ik = 0) ~ l/(g0o). Since (fio ~ 1/g the amplitude of the quantum fluctuations 
are always 0(1). 

It is clear that as we vary g the energy density of the system which is of order g 2 <pQ 
varies as well. In order to see the parametric behavior of the relaxation time on g at 
fixed energy density we instead take 0o ~ l/<7 1 • in this case we find that 

r rclax ~ / ~ -^7T (29) 

9y/V*(k = 0) 9 m 

which is consistent with the 1/g 2 / 3 behavior extracted from the 3D simulation of |20j . 
In Fig. [6] we show the number density defined by 

n k = <0|ota k |0> = ~ ^ k |0 k | 2 + - i (30) 

at various times along the evolution. The initial condition is such that the number 
density is zero, n k = 0, except for the zero-mode which is highly occupied. As the 
system evolves one sees the appearance of peaks. Whether these peaks correspond to 
resonance modes is not clear. Even though we know the location of the resonance band 
at t ~ the effective mass of the background field changes with time thereby changing 
the location of the resonance band with time. Of course, larger lattice simulations will be 
needed to reinforce these statements. At late times, when the system has fully relaxed, 
the number density is smooth with a power-law fall off. 

In Fig. 0we show the final number density. In this case we have plotted the spectra 
at the discrete values allowed by our grid. The solid curve is a fit to n k ~ with 
s = 1.45. It clearly does not fall as l/w k as one would expect from classical thermal 
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Figure 6. Evolution of the number density at various times. The curves at t = 160 
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Figure 7. Final number density versus momentum. The solid curve is the best fit to 
~ wr 11 with s = 1.45. 



equilibrium. Interesting further work might study the late time behavior of the particle 
number to see if it scales according to Kolmogorov turbulence J3TJ [39] . 



6. Conclusions 



In conclusion, it is apparent that quantum fluctuations modify the evolution of a classical 
scalar theory to the point where it evolves in accordance with ideal hydrodynamics. We 
have attempted to understand this behavior by using modified spectra of fluctuations 
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and by studying the particle content of the theory. We observe that there is a two-stage 
relaxation process; the rapid early time relaxation is somehow controlled by modes of 
intermediate momentum (at and above the resonance band) followed by a longer more 
gradual relaxation which is controlled by the lower momentum (near zero) modes. The 
methods used for the scalar field can presumably be extended to the case of gauge 
theories. 
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